Chiral Symmetry Restoration from a Boundary 



B. C. TiburzE 

Department of Physics, The City College of New York, New York, NY 10031, USA 
Graduate School and University Center, The City University of New York, New York, NY 10016, USA 
RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973, USA and 



o 



> 

in 

O 

m 



Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106, USA 

(Dated: February 28, 2013) 

The boundary of a manifold can alter the phase of a theory in the bulk. We explore the possibility 
of a boundary-induced phase transition for the chiral symmetry of QCD. In particular, we investigate 
the consequences of imposing homogeneous Dirichlet boundary conditions on the quark fields. Such 
boundary conditions are employed on occasion in lattice gauge theory computations, for example, 
when including external electromagnetic fields, or when computing quark propagators with a reduced 
temporal extent. Homogeneous Dirichlet boundary conditions force the chiral condensate to vanish 
at the boundary, and thereby obstruct the spontaneous breaking of chiral symmetry in the bulk. 
As the restoration of chiral symmetry due to a boundary is a non-perturbative phenomenon, we 
utilize the sigma model to exemplify the issues. Using this model, we find that chiral symmetry is 
completely restored if the length of the compact direction is less than 2.0 f m. For lengths greater 
than about 4 f m, an approximately uniform chiral condensate forms centered about the midpoint of 
the compact direction. While the volume- averaged condensate approaches the infinite volume value 
as the compact direction becomes long, the finite-size corrections are shown to be power law rather 
than exponential. 

PACS numbers: 12.39.Fe, 12.38.Gc 



I. INTRODUCTION 

One of the hallmark non-perturbative features of QCD 
is spontaneous breaking of chiral symmetry. While quark 
masses explicitly break the chiral symmetry of the ac- 
tion, the lightest quarks (up and down) have masses that 
can be treated as a perturbation about the symmetric 
SU(2) L <g) SU(2) R chiral limit. The formation of a chi- 
ral condensate by the QCD vacuum in the chiral limit, 
namely < ipx/j >^ 0, spontaneously breaks the chiral 
symmetry down to the vector subgroup, SU(2)y. This 
symmetry breaking pattern along with the explicit break- 
ing due to the quark masses gives an explanation of the 
lightness of the isotriplet of pseudoscalar pions because 
they must be the emergent Goldstone bosons. 

Lattice gauge theory provides a first principles method 
for solving QCD numerically on finite Euclidean space- 
time lattices. Strictly speaking, spontaneous symmetry 
breaking cannot occur in a finite volume [lj. In prac- 
tice, the formation of a chiral condensate on periodic 
lattices is determined by the size of the pion Compton 
wavelength compared to the lattice size [2]. When the 
pion Compton wavelength is small compared to the lat- 
tice size, l/m n -C L, the finite volume effect on the con- 
densate is exponentially small and can be computed in 
chiral perturbation theory [3|. On the other hand, when 



the pion Compton wavelength is large compared to the 
lattice size, l/m^ 3> L, the zero-momentum mode of the 
coset of Goldstone fields becomes strongly coupled 
Taking the chiral limit at finite volume puts one in the 
latter regime, which is the regime in which chiral sym- 
metry is restored. 

In this work, we explore a different restoration of chiral 
symmetry!]] We consider the fate of chiral symmetry on 
a Euclidean manifold with three infinite directions, and 
one compact direction. The compact direction is sub- 
ject to homogeneous Dirichlet boundary conditions. Such 
boundary conditions are employed on occasion in lattice 
gauge theory computationso One example occurs in the 
study of hadron properties in external electromagnetic 
fields. Naive implementation of uniform external fields 
via linearly rising four- vector potentials leads to field gra- 



' btiburzi@ccny.cuny.edu 



Throughout we refer exclusively to the non-singlet chiral sym- 
metry, SU(2)i J ® SU(2)r. For the singlet case, it is convenient 
to phrase the discussion in terms of the axial symmetry of the 
QCD action, U which is anomalous in infinite volume. On a 
manifold with a boundary, the pressence of the axial anomaly is 
subtle 5]. With a homogeneous Dirichlet boundary, it is known 
that the integral of the singlet axial current's divergence van- 
ishes [0| . We leave the investigation of the chiral anomaly to 
future work. 

Inhomogcncous Dirichlet boundary conditions are imposed in 
the Schrodinger functional representation of QCD [jj . These 
boundary conditions too can affect spontaneous chiral symmetry 
breaking as will be made clear in a footnote below. 
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dients at the lattice boundary, and homogeneous Dirich- 
let boundary conditions have been sought to mitigate the 
effect of these gradients on the quark fields pMLlj. Ad- 
ditionally temporal Dirichlet boundary conditions have 
been imposed on lattices to compute quark propagators 
with a reduced temporal extent (so-called chopping of 
lattices) [IM3I3 

The possibility that imposing boundary conditions can 
affect the phase of a theory in the bulk is a known phe- 
nomenon in condensed matter physics, for an overview of 
boundary critical phenomena, see [181 ]. There has been 
some very recent work using the Gross-Neveu model to 
investigate the relation between chiral symmetry break- 
ing and the Casimir force through a direct computation 
of the boundary separation and tem per ature dependence 
of the model's vacuum energy fl9l. |2C|. To our knowl- 
edge, the effect of homogeneous Dirichlet boundary con- 
ditions on the chiral condensate in QCD has otherwise 
not been explored. We employ a different approach to 
demonstrate that homogeneous Dirichlet boundary con- 
ditions can lead to the restoration of chiral symmetry. 
Our consideration is in the absence of any external elec- 
tromagnetic fields. Unlike on periodic lattices, we find 
that the restoration of chiral symmetry in the case of 
Dirichlet boundary conditions is not controlled by the 
size of the pion Compton wavelength. Instead, it is the 
Compton wavelength of the sigma meson which strongly 
affects the restoration of chiral symmetry in the pres- 
ence of a Dirichlet boundary. To investigate these effects, 
we utilize the sigma model and constrain the model pa- 
rameters from phenomenology. We find substantial ef- 
fects on the chiral condensate for lattice sizes less than 
about 4fm. As our consideration lies outside of chiral 
perturbation theory, this estimate is necessarily model 
dependent. The sigma model, however, has the minimal 
features necessary to address the effect, whereas chiral 
perturbation theory by contrast does not. An important 
finding within this model is that finite-size effects on the 
volume-averaged condensate only vanish as a power law 
rather than as an exponential. 



In both of these scenarios, only the valence quarks are subjected 
to Dirichlet boundary conditions, and one should formulate the 
problem in terms of graded symmetries [l5, 16] . For a partially 
quenched theory utilizing the same boundary conditions in the 
valence and sea sectors, chiral symmetry breaking takes the form 
SU(4\2) L <g> SU(4\2) R -> ST/(4]2)v. With different boundary 
conditions, however, there is no symmetry relating the valence 
and sea sectors, and the chiral symmetry of the action must have 
the form SU(2\2) L ® SU(2) L ® SU(2\2) R <g> SU(2) R , which is the 
same symmetry as in mixed-action lattice QCD calculations [17| . 
This reduced symmetry can break either to SU (2\2)y (g> SU {2)y, 
by the formation of valence and sea quark chiral condensates, 
or to 1 ® SU(2)v, by the formation of only a sea quark chi- 
ral condensate. Our consideration above is for the valence sub- 
group of the graded chiral symmetry, i.e. SU(2)l <%> SU(2) R C 
SU{2\2)l ®SU(2\2) r , and the mesons of the sigma model should 
be viewed as built only from valence quarks. One can go fur- 
ther and generalize the sigma model to appropriately reflect the 
graded symmetries; but, this is beyond the scope of our work. 



Our presentation has the following organization. We 
first describe the simple sigma model that we employ 
in Sec. [Til The model is then considered in Sec. IIIII with 
homogeneous Dirichlet boundary conditions in one of the 
space-time directions. We obtain an expression for the lo- 
cal chiral condensate by solving the analytical mechanics 
that determines the vacuum. The volume-averaged con- 
densate is shown to approach the infinite volume value 
up to power-law corrections. Finally in Sec. IIV1 a brief 
summary concludes this work. 

II. SIGMA MODEL 

We begin by describing the simple sigma model we 
employ. The ingredients are chosen so that the model 
shares the same pattern of symmetry breaking as QCD 
with two light flavors. Our consideration in this section 
is restricted to infinite volume. The model consists of 
an isoscalar scalar field, S, and an isovector pseudoscalar 
field, P. These fields appear in the Euclidean space ac- 
tion density having the form 

C E = \d^S d^S + ±d„P ■ d^P - 

+k(s 2 + P 2 -v 2 ^ . (1) 

There are four parameters of this model, m, v, A, and 
A. The parameter m is analogous to the quark mass. 
When the quark mass parameter is set to zero, the ac- 
tion density in Eq. (JD has an 5*0(4) ^ SU(2) ® SU(2) 
symmetry. The quartic interaction term, however, leads 
to spontaneous breaking of the 5*0(4) symmetry down 
to 5*0(3) = SU(2), by the formation of a condensate, 
namely S 2 + Pq = v 2 . We append a subscript "0" 
throughout to denote the vacuum expectation value of 
a field. Here the vacuum expectation values are uniform 
in space-time, e.g. Sq(x) = Sq. 

In discussing the sigma model, it is convenient to use a 
polar decomposition of the fields. To this end, we write 

S + iP ■ f = E{7, (2) 

where r are the usual isospin matrices, S is a real-valued 
field, and U is a unitary field. The latter encompasses 
the Goldstone bosons if in the form U = exp(i7r • t/F). 
In terms of the polar decomposition, the sigma model 
action density has the form 

C E = -tr [d M £d M £ + tfdnUdpUi] 

~ix[mE(U + U*)] +A(£ 2 -w 2 ) 2 . (3) 

Traces appearing above are taken over isospin. The vac- 
uum state of the theory is determined by minimizing the 
action density. The action has a minimum for the values Q 



4 In the chiral limit, one has two possibilities for the vacuum 
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Uq = 1 and So = v. Expanding about these vacuum ex- 
pectation values, U = 1 + + • • • , and £ = w + er, we 
have to quadratic order in the fields 

C = ^d^ad^a + \^n 2 a a 2 + ^d^n ■ + \^m\ n • tt, (4) 

with the meson masses identified as m 2 = 8AF 2 , and 
m 2 — Xm/F 2 . Notice we must have the equality F = v 
in order for the pion kinetic term to have the proper 
canonical normalization. 

The chiral condensate, which we denote by < tjnp >, 
can be found by differentiating the vacuum energy den- 
sity with respect to the explicit symmetry breaking pa- 
rameter, which is the quark mass to. The value of the 
chiral condensate is simply 

< rfip >= -A. (5) 

This identification enables us to rewrite the pion mass 
in the form F 2 m 2 — \ < tpip > | to, which is the Gell- 
Mann-Oakes-Renner relation. In the absence of explicit 
symmetry breaking, i.e. the chiral limit when m = 0, the 
pions are massless Goldstone modes. The sigma meson, 
on the other hand, remains massive, with its mass scale 
set by the pion decay constant F. In the chiral limit, 
the sigma meson can thus be integrated out to pro duce a 
low-energy theory of just the Goldstone modes [2 1] . This 
theory is chiral perturbation theory and is a model inde- 
pendent description of low-energy QCD [22]. To study 
the effects of a Dirichlct boundary on the chiral conden- 
sate, however, we must retain the sigma degree of free- 
dom. Without the sigma, the chiral condensate would be 
frozen into a uniform value, and this conflicts with the 
quark boundary conditions. 

Finally we can use phenomenology to fix the model pa- 
rameters. With our normalization, the vacuum expecta- 
tion value v of the sigma field is identical to the parame- 
ter F, which can be identified as the pion decay constant 
once an external axial- vector field has been coupled to 
the theory. Technically F is the chiral limit value of the 
pion decay constant, and has the approximate numeri- 
cal value F = 93 MeV. The mass of the sigma resonance 
has been determined from a detailed analysis of high- 
precision tttt scattering data (23[. From the extracted 
mass of the sigma, m a = 440 MeV, we find the parame- 
ter A has the value A = 2.8. In our model, the sigma 
meson is a particle not a resonance, and consequently 
the usual caveats about meson models apply. The chiral 
condensate can be used to fix the value of the parameter 
A, which itself is the chiral limit value of the condensate 
up to sign. This quantity, however, is QCD renormaliza- 
tion scale and scheme dependent, and we simply choose 
to quote results in terms of ratios to the infinite volume 



expectation value of the sigma field, So = ±f . Inclusion of a 
small quark mass term, and the subsequent chiral limit, m — > 0, 
will lead to the selection of Eq = +v. 



condensate, for which the scale and scheme dependence 
exactly cancels 



III. DIRICHLET BOUNDARY 

Above we describe the symmetry breaking in the sigma 
model in infinite volume. To minimize the action as a 
functional of the meson fields, we eliminated the kinetic 
terms by restricting our attention to uniform field config- 
urations. These uniform field configurations are energet- 
ically preferred because the kinetic energy contributes 
positively to the action density. With a homogeneous 
Dirichlet boundary condition, however, there are no non- 
trivial field configurations that are uniform. On the other 
hand, the potential energy terms in the action drive the 
£ field to a non-vanishing vacuum expectation value. In 
general, the Dirichlet boundary now allows for compe- 
tition between the kinetic and potential terms, and the 
symmetry breaking pattern of the model must be scruti- 
nized. 

For simplicity, we choose to impose a homogeneous 
Dirichlet boundary condition in just one of the space-time 
directions. Because we work in Euclidean space, we need 
not specify whether it is a spatial or temporal direction 
that has a boundary. The remaining three directions 
are kept infinite. We thus consider the quark field ip(x) 
to satisfy the boundary conditions ip(x = 0) = 0, and 
ip(x = L) = 0, where L is the length of the x-direction0 



5 It is efficacious to spell out a procedure to investigate the effects 
of Dirichlet boundary conditions on the chiral condensate using 
the lattice as a regulator. In this case, the condensate is both 
multiplicatively and additively renormalized at finite values of 
the quark mass. We can define a subtracted chiral condensate for 
the case of Dirichlet boundary conditions in a way that is anal- 
ogous to the temperature dependent case, see, for example, |24| . 
Consider two values of the quark mass, mi, and m,2 for a given 
lattice action. Using this ultraviolet regulator, one computes four 
quantities: the chiral condensate with periodic boundary condi- 
tions for each quark mass, < tf>tp > mi and < ipip >m 2 > an d the 
local chiral condensate with Dirichlet boundary conditions for 
each quark mass, < iptp(x) > mi and < iptp(x) >m 2 - The sub- 
tracted condensates, < iptp > su b=< ipip >m! < V"/' >m 2 
and < ifnpix) > su b = < ipipix) > mi < i>ip(x) > m2 , are free 
of power-law divergences, which are proportional to ma~ 2 . In 
the latter case, the subtraction cancels a quark-mass independent 
function that depends on x. The ratio of subtracted condensates, 
< ipip(x) > su b / < ipip >subi can then be used to study the be- 
havior of the chiral condensate in the infrared. Because the chiral 
condensate calculated with periodic boundary conditions suffers 
from finite-size effects, one should first take the infinite volume 
limit of < ipip > su b in order to isolate the finite-size effects aris- 
ing solely from Dirichlet boundary conditions when forming the 
ratio. Finally the local chiral condensates can be averaged over 
the compact direction to investigate the behavior of the volume- 
averaged condensate. Such quantities are more likely to be com- 
puted using the lattice in order to reduce statistical noise. 

6 In the Schodinger functional representation of QCD Q, one im- 
poses inhomogeneous Dirichlet boundary conditions in the tem- 
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The other space-time directions are treated implicitly in 
our notation. 

The imposition of quark boundary conditions trans- 
lates into boundary conditions on the E and 7r fields 
of the sigma model. These meson fields have the same 
quantum numbers of certain quark-level interpolating op- 
erators. In particular, we have £(x) ~ ip{x)^{x) and 
7?(x) ~ ip(x)j5Tip(x). As a result of these identifications, 
we have the boundary conditions 



£(x = 0) = £(x = L) = 0, 
U(x = 0) = U(x = L) = 1. 



(6) 



These boundary conditions can be directly obtained if 
one derives the sigma model from a model with quarks, 
namely from the Nambu-Jona-Lasinio model [27[ . 

To find the vacuum of the theory with Dirichlet bound- 
ary, we minimize the action density. For the directions of 
infinite extent, the fields are frozen into their correspond- 
ing zero momentum modes. On account of Eq. (jf)]), the 
unitary field U can take on a uniform value Uq = 1 for 
all space-time. The quark mass term will consequently 
be minimized as a function of Uq for Uq = 1 provided 
the vacuum value of the field £ is positive. As the quark 
mass dependent case is more complicated, we will work 
in the chiral limit, m = 0, for which the vacuum mini- 
mization requires Uo = 1. It is important to note that 
adding a small quark mass term to the action leads to an 
energetic preference for positive values of £. 

The vacuum value of the field £ depends on the in- 
coordinate, namely £ = £o(x) + • • • . In the chiral limit, 
this vacuum expectation value will minimize the action 
density 



M£nl = I 



dx, (7) 



which is a functional of £o(x). The functional minimiza- 
tion, SS/6'Sq(x) = 0, can be achieved through solving the 



poral direction having the form, 'P+V( :r 4 = 0) = i/j(x4 = 0)'P— = 
0, and V-ipixA =T)= ^(x 4 = T)V+ = 0, where V± = §(1±7 4 ) 
are parity projection matrices, and the remaining components of 
the fermion field are non-vanishing at the boundaries. From a 
renormalization group argument, such boundary conditions are 
natural in the continuum limit 12511 . The chiral condensate can 
be written trivially as the sum, < ipip(x) >=< ipV+ip(x) > 
+ < tp'P—ip{x) >, which consequently vanishes at the boundaries 
X4 = 0, and T. The Schrodinger functional is largely employed as 
a renormalization scheme; and, as a massless scheme, the mass- 
less fermions cannot scatter off the chiral condensate. Restora- 
tion of chiral symmetry due to the boundary is then irrelevant 
for determining the renormalization of operators in this scheme. 
The same is not true for the calculation of hadron properties at 
non-vanishing quark masses within the Schrodinger functional 
formulation. While the utilization of small temporal lattices be- 
comes possible to compute hadron properties, see, for example, 
[26l | , one must be careful to account for possible non-perturbative 
finite-size effects resulting from the temporal dependence of the 
chiral condensate. 



Euler-Lagrange equation. Once this solution is known, 
the value of the chiral condensate follows immediately 



< tj}ip(x) >= ~J^{x). 



(8) 



Accordingly the chiral condensate will vanish at the 
boundary, but could develop a non-zero value in the bulk. 
It remains to solve the Euler-Lagrange equation to deter- 
mine £o(x). 

The problem of minimizing the action density is equiv- 
alent to a problem in analytical mechanics. The mechan- 
ical analogue of energy is given by 



£ 



d£o 

dx 



(9) 



and is a constant of the motion. Because of the Dirichlet 
boundary at x = and x — L, any non-trivial solution for 
£o(x) must have at least one turning point. The turning 
points Xj are determined by the condition d£o /dx | _ = 
0. Labeling the corresponding values of the field at the 
turning points, Eq = Eo(xj ), we see that the latter are 
given by 



^(3) 



-,(2) 
.(4) 



= +vy/l + t, (10) 



with £ = sj-^s. In order for there to be any turning 
points at all, it must be the case that £ < 0. For the 
first two turning points to exist, one additionally requires 

£<i. 

The concavity of the function Eo(x) follows from the 
analogue of Newton's second law. The force is given by 
F = -dV/dT, , with V = -A(Eg - v 2 ) 2 , and implies 

that the values of the field at the turning points, Eq , 
are local minima for j even, and local maxima for j odd. 
As a result, we need not consider the turning points E 

and £ 4 ^ in determining the solution Eo(x). For exam- 
ple, consider a solution which rises from £0 = at x = 0. 
For this solution to turn over, we need a positive turn- 
ing point corresponding to a maximum. The only pos- 
sibility is the value Eq 1 ^. The solution Eo(x) could then 
decrease from the maximum down to a minimum before 
rising again, but the only possible turning point that cor- 
responds to a minimum with value less than Eq 1 ^ is £q 2 ^- 

The solution £o(x) is necessarily bounded by £q from 

(2) 

above and £ from below. There are an infinite num- 
ber of solutions which minimize the action. They are 
characterized by the number of oscillations between the 
extrema Eg 1 ' 1 and E Q 2 ^. Because the latter is negative, the 
inclusion of a small quark mass term will lead to an en- 
ergetic disadvantage for all solutions having any turning 
points for which £o(x) attains the value E . 

The solution Eo(x) we seek can thus be characterized 
as monotonically increasing from zero up to the maxi- 
mum value Eq 1 ^, and then monotonically decreasing back 
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L [fm] 



FIG. 1. The solution £ as a function of the finite extent L 
of the ^-direction. Values of L for which £ > 1 lead to the 
vacuum expectation value Eo(x) = 0, and hence correspond 
to a complete restoration of chiral symmetry. 




x/L 

FIG. 2. Ratio of the chiral condensate with homogeneous 
Dirichlet boundary conditions to the infinite volume chiral 
condensate plotted as a function of the ^-coordinate scaled 
by L. For L < 2.0 fm, the condensate vanishes everywhere, 
< ^ip(x) >= o. 



down to zero. The motion of So is symmetric about the 
turning point; consequently, we must have x\ = L/2. In- 
tegrating the equation of motion from the boundary to 
the turning point (or vice versa), we arrive at the equa- 
tion 



1-g 



■6 



11) 



where K(to) is the complete elliptic integral of the first 
kind. This is a special case of the incomplete elliptic 
integral of the first kind, F(cf> |m), defined by 



> |m) = 



dO 



o yT 



(12) 



to sin 



namely K(to) 
abbreviation 



F(-||m). Above, we have employed the 



I = vLV2A = ^m a L. 



(13) 



The relation expressed in Eq. (jllj) implicitly defines the 
analogue of the mechanical energy £ as a function of 
the sigma model parameters and the extent of the x- 
direction, i.e. £ = £(v, A,L), In practice, it is simpler to 
work with the dimensionless variable £ ~ £(v,A,L). In 
Fig. [TJ we plot the value of £ that satisfies Eq. (ITT1) as a 
function of the extent of the compact direction, L. For 
L < - ^= = 2.0 f m, the solution requires £ > 1 for which 

the turning point does not exist, and consequently 
chiral symmetry is completely restored. 

With the value of £ determined from Eq. (jlip , we can 
implicitly specify the solution Eo(x) by integrating the 
equation of motion to an arbitrary point x. We find the 
solution must satisfy 



sin 



.(i) 



l-£ 



1 + C 



L ■ 



for 
for 



< x < § 
-§ < x < L 



(14) 



with F(i/)|to) the incomplete elliptic integral defined 
above. Notice the full solution agrees with Eq. ([TTj) at 
the turning point. 

From this solution, we can determine the chiral con- 
densate as a function of the x-coordinate. It is simplest 
to consider the ratio of the condensate with Dirichlet 
boundary conditions with that of the condensate in infi- 
nite volume, namely < ipip(x) > / < ipip >■ This ratio 
is plotted in Fig. [5] for several values of the extent L. 
For L < 2.0 fm, the chiral condensate identically van- 
ishes for all x. Above this size, a non-vanishing con- 
densate forms, however, the condensate is significantly 
altered from its infinite volume value for sizes less that 
about 4 f m. When L = 3.5 f m, for example, the conden- 
sate is reduced by more than 25% over half of the length 
of the x-direction. The condensate, moreover, reaches 
a maximum value that is 15% smaller than the infinite 
volume value. For sizes greater than about 4 f m, an ap- 
proximately uniform chiral condensate is formed. The 
support of the condensate, however, is over roughly 50% 
of the extent of the x-direction, and is naturally centered 
about the midpoint x — L/2. 

While the chiral condensate must vanish at the bound- 
ary for all L, chiral symmetry is spontaneously broken in 
the bulk of the lattice in the limit of large L. The maxi- 
mum value of the chiral condensate approaches the infi- 
nite volume value, because the sigma field at the turning 
point has the behavior, Eq — > v as L goes to infin- 
ity. This behavior arises due to £ approaching zero, see 
Fig- ID From Eq. ([TTj) . we can analytically derive the large 
L limit of the maximum value of the condensate. Us- 
ing the expansion of the complete elliptic integral about 
unity, we find the midpoint value is given by 



< ijnp(L/2) >=< i/itp > 



1 — 4 exp 



which shows that the asymptotic behavior is controlled 
by the Compton wavelength of the sigma meson. 
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FIG. 3. Ratio of the volume-averaged condensate < ipij) > to 
the infinite volume condensate < ifjtp > plotted as a function 
of the finite extent L. The dotted curve shows the asymptotic 
formula given in Eq. (|17[) . 

The value of the chiral condensate at the midpoint does 
not completely characterize the fate of chiral symmetry 
because there is also the issue of support over the bulk 
of the lattice. In addressing the extent to which chi- 
ral symmetry is restored, one can also study the chiral 
condensate averaged over the compact direction. This 
average is simply defined by 



< ipifj > = — dx < ipip( x ) > ■ ( i6 ) 
L Jo 

In the limit of an asymptotically large extent L, the 
volume-averaged chiral condensate tends to the infinite 
volume value, however, the approach to asymptopia is 
slow. In the asymptotic limit, we expand the second 
argument of the incomplete elliptic integral in Eq. (fT4|) 
about unity. The asymptotic condensate can be deter- 
mined in closed form, and averaged over the ^-direction 
to produce 

<l/o=<#> (l-^ + --- ]. (17) 

This power-law behavior is confirmed in Fig. [31 where 
we compare the volume-averaged condensate determined 
from the full solution, Eq. (JUJ), to the asymptotic form 
given in Eq. (|17[) . While the asymptotic form works very 
well for L > 3 f m, the volume-averaged condensate only 
slowly approaches the infinite volume value. 

IV. SUMMARY 

Above we investigate the effect of homogeneous Dirich- 
let boundary conditions on the chiral condensate. This 
effect cannot be ascertained within chiral perturbation 
theory, because the chiral condensate is determined by 
the expression 

<4>ip{x) >= -~ < U(x) +U^(x) > +■■■ . (18) 



The right-hand side of this relation does not vanish at 
the boundary in contradiction with the quark bound- 
ary conditions satisfied on the left-hand side. A correct 
treatment of the chiral condensate in the presence of a 
Dirichlet boundary necessitates including the dynamics 
of the sigma meson. For this reason, we employ the sigma 
model which shares the same symmetry breaking pattern 
as QCD with two light quark flavors. Model parameters 
are chosen to reproduce phenomenology, and we subject 
the model to homogeneous Dirichlet boundary conditions 
in one of the space-time directions. The chiral conden- 
sate is determined by minimizing the action to find the 
vacuum configuration, and this is achieved in the sigma 
model using well-known functions. 

Chiral symmetry is shown to be restored in the sigma 
model for sizes less than 2.0 fm. For sizes greater than 
about 4 f m, an approximately uniform condensate forms 
centered about the midpoint of the compact direction 
and having support over roughly half of the compact di- 
rection. The volume-averaged condensate suffers from 
finite-size effects that are power law rather than expo- 
nential. As our estimation is necessarily model depen- 
dent, it would be useful to compare the behavior of the 
chiral condensate in quark models with spontaneous chi- 
ral symmetry breaking, and other meson models with 
resonances. 

Nonetheless, we believe one should be cautious in in- 
terpreting results from lattice computations employing 
homogeneous Dirichlet boundary conditions on lattices 
less than about 4fm in extent. As our estimate, more- 
over, comes from treating only one space-time direction 
as compact, it would be interesting to consider the case 
of two directions subject to Dirichlet boundaries. Such 
a scenario is encountered in some lattice computations 
with external electromagnetic fields (2|| [29| . This setup 
would additionally allow for an exploration of the Casimir 
effect. The three-dimensional case is also of interest to 
confront phenomenology of the bag model [30| l3l| , and 
the proposal of in-hadron condensates [lU, [33| . Finally 
to establish the credibility of lattice computations with 
homogeneous Dirichlet boundaries, one requires a lattice 
computation of the chiral condensate, either locally or 
volume averaged, which will ultimately reveal the extent 
to which chiral symmetry is restored in the presence of a 
boundary. 
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